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Abstract 

^ We present next-to-leading order QCD corrections to production of two W bosons in 

■ hadronic collisions in the extra dimension ADD model. Invariant mass and rapidity 

distributions are obtained to order as in QCD by taking into account all the parton level 
subprocesses. The computation is organized using the monte carlo based method of phase 
space slicing. We estimate the impact of the QCD corrections on various observables and 
find that they are significant. We present some results for a 10 TeV LHC but most of the 
results presented here are for 14 TeV LHC. We also show the reduction in factorization 
scale uncertainty when 0{as) effects are included. 
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1 Introduction 



The fact that electroweak symmetry breaking scale of the Standard Model (SM) cannot 
be made stable against the quantum corrections (hierarchy problem) within the SM 
indicates to the possibility of new physics at TeV scale. The Large Hadron Collider 
(LHC) which will operate at an enormous center of mass energy (v^ = lATeV) offers 
to shed light on the existence of new physics. The most popular new physics models are 
based on the ideas of supersymmetry and extra spatial dimensions. Proposals to address 
the hierarchy problem using extra dimensions were introduced in [HE]- In this paper 
we will consider the model by Arkani-Hamed, Dimopoulos and Dvali (ADD) [Ij. In this 
model all SM fields are confined to a (3+1) dimensional manifold and the extra d spatial 
dimensions are compactified, with same radius of compactification R, on a (i-torii. The 
effect of extra dimensions appears as Kaluza Klein (KK) gravitons on the 3-brane which 
couple to SM fields through energy momentum tensor with a strength k, which is related 
to the volume of the extra dimensions, and the fundamental scale Mg in 4 + ci dimensions 
by [3] 

K^R'^ = 8n{AnY/^T{d/2)M-^'^+^\ (1) 

Although the coupling k is Mpi suppressed, the fact that there are large number of KK 
modes that couple to the SM fields makes the cumulative effect significant and leads to 
observable effects. One extra dimension ie. d = 1 is ruled out |1] and d = 2 is severely 
constrained so we will consider in this paper d = 3 and above. There are two ways to 
probe such effects at colliders, either through gravitons emission or by virtual graviton 
exchange. In this paper we will consider only the effects of virtual spin-2 KK states. 

The precise measurement of hadronic production of gauge boson pairs is one of the 
important endeavors at the LHC both in the context of SM and new physics studies. 
Studies in other channels have been reported in [5] in extra dimension models. In this 
paper we will consider production of W pair at the LHC. Owing to its importance, 
its study has attracted a lot of attention in the literature. Many studies have been 
carried out for its production in the SM; a study in the context of anomalous triple 
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gauge boson vertices was carried out in [6]. Leading order (LO) studies in the SM can 
be found in [7j. As is well known the LO results are highly sensitive to the arbitrary 
renormalization and factorization scales. At this order the factorization scale enters 
solely through the parton distribution functions as the parton level cross-section, at this 
order, does not depend on fip- As we include higher order terms of the perturbation 
series the dependence will reduce and an all order result will be completely independent 
of these arbitrary scales. In addition the NLO results are usually significantly enhanced 
as compared to the LO results. It is thus important to carry out NLO calculation to 
reduce sensitivity to these scales. Because of its importance, W^W~ production has 
been studied to next-to- leading-order (NLO) accuracy in the SM [8l[9]. It has also 
been studied via gluon fusion through a quark box loop or triangle quark loop with 
7 or Z boson exchange [lO] and at one and two loop levels in high energy limit in 
SM [llj. The significance of NLO computations in the extra dimension models for Drell- 
Yan [12j, diphoton [13], ZZ p[4], graviton+photon [15j, graviton+jet [16j production has 
already been demonstrated. These studies show that not only the predictions at NLO 
are enhanced but are also less sensitive to the factorization scale. With this in mind we 
carry out a complete NLO calculation of hadronic W^W~ production in ADD model in 
this paper. W^W~ production in Randall Sundrum model is presented in our work jlTj . 

We organize the paper in the following sections as follows. In section 2 we give the 
details of NLO computation for the W^W~ production. Here we give the matrix element 
squares for all the 2 — 2 subprocesses both at leading order and at loop level at order 
as (= (?f/167r^). All the Feynman diagrams are collected in the appendix. The virtual 
corrections contain soft and collinear singularities which appear as simple and double 
poles in e as we use dimensional regularization with n = 4 + e. The 2—7-3 real emission 
matrix elements which were also calculated analytically are not presented here (to save 
space as the expressions are lengthy) and can be obtained on request. These compu- 
tations were done using the symbolic manipulation program FORM [18|. We intend to 
use monte carlo methods for obtaining kinematical distributions. These methods prove 
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very useful if experimental cuts need to be imposed on the final state particles and they 
avoid the need for repeating calculations for obtaining different distributions as would 
be required by completely analytical methods. In this paper we employ the method 
of two cutoff phase space slicing [19j to carry out our calculation. This method gives 
numerically stable results as has been demonstrated in [19j and also in our previous 
works [inilTl]. An alternative to this is dipole subtraction method which is also 
widely used in higher order QCD calculations. We describe, in brief, phase space slicing 
method and how the soft and coUinear singularities that appear at virtual level and real 
emission level are treated. In section 3 we present some checks on our code and then 
present some useful kinematical distributions. 

2 Next-to-leading-order Computation 

The hadronic production of W bosons at NLO has three pieces of computation. A LO 
piece which is a 2 —t- 2 parton level process {qq — j- W~^W~ and gg — > W'^W~); second 
is the 2 — 2 order piece which originates from loop corrections to the LO piece; the 
third and final part originates from real emission processes where in addition to two W 
bosons, a parton is also emitted in the final state {qq — )■ W~^W~g, qg — )■ W~^W~q, and 
gg — )■ W^W^ g ). Let us take up these three pieces in turn. 

2.1 Born Contribution 

The charged vector-boson production in the leading partonic scattering processes corre- 
sponds to 

<Pi) + h{p2) ^ W^{p^) + W-{p^). (2) 

where pi and p2 are the momenta of initial state partons while ^3 and p4 are those of final 
state vector bosons. The W boson pair can couple to KK gravitons, so it is possible to 
produce them through virtual graviton exchange at the leading order [21]. In SM this 
proceeds via t channel (or u channel) quark anti-quark annihilation along with s-channel 
7 and Z boson exchange shown in Fig. [TOl The coupling of fermions to W and Z bosons 
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are respectively 



- (1 - 7^) , -^'^1' {C^ - Cal') . (3) 



where 7>v and Tz read 

V2 sm 6^11/ sm 6w cos fc^vK 

The coefficients and Ca are given by 

a = - 2 Q/ sin2 ^i^, Ca = Tl (5) 

The Qf and T/ denote the electric charge and the third component of the weak isospin of 
the fermion / respectively, 6w is the weak mixing angle and is the mass of Z boson. 
The Z boson propagator is given as 



s — ml + i TzTUz 

We have used unitary gauge in the electroweak sector ie. ^ = oo, this simplifies the 
calculation as both the goldstone boson and ghosts in the electroweak sector disappear. 

The 75 matrices that appear in the intermediate stages of the computation require 
special care as they are not defined in arbitrary dimensions. We have used naive anti- 
commutation relations between 75 and other gamma matrices in n dimensions and the 
resulting traces are then computed in n dimensions as they are free of 75. Alternatively, 
one can use other method namely HVBM-scheme which was proposed in [22] and gen- 
eralized in [23]. In this approach. Gamma matrices and momenta in the loop and final 
state phase space integrals are split into a 4 and an n — 4 dimensional part. The 75 
anti-commutes in 4 dimensions and commutes in n — 4 dimensions with rest of the 7 
matrices. 

We give below the matrix element squares summed (averaged) over the final (initial) 
state spins, colors and polarizations notated by an overline in n = 4 -|- e dimensions . 
s, t and u are the usual Mandelstam invariants. We will denote by sm, gr and int the 
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contributions from SM, gravity, and interference of SM with gravity respectively. The 
SM at LO gives order contribution to the cross sections: 

qq,sm 



|M(o)|2 



N 



{AlBf + AlBl + AlBl) 



(7) 



where 



Al 



'W ' 



a 

-^ + eqCyTz cot 6w 



[s — mz 



s [{s — miy + r^m^] 
1 



AU 

^3 



+ (C„ + a) 7^ Tz coiOw 



2 [(. - m2)2 + r>2] 



(8) 



Al are in their essence combinations of EW couphngs and propagator factors and Bl are 
the functions of the kinematic invariants. N denotes the number of colors. In case of 
up type quark initiated processes, -B" originates from purely t-channel, B2 from purely 
s-channel while B^ from the interference of t and s-channel diagrams. These functions 
are given by 

1 



B't 



- [m^ (-2 + n) ] + t^u - 2m^ (-2 + n) (t + u) 



+mH {{-9 + An) t + {-2 + nf u} 



(9) 



2m4 



- 20m® (-2 + n) + 8m^ (-2 + n) {t + u) + tu {t + u) 



-2m^ (-2 + n){t + uf + (-9 + 4n) t'^ + 2 (-13 + 6n) 
+ (-9 + 4n) 



tu 



(10) 
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2mH 



- lOm^ (-2 + n)+ Am^ (-2 + n) {t + u) + (t + u) 



-2m^ (-2 + n)t{t + uf + mH{{-9 + 4n) t + (-13 + 6n) u} 



'111 



Here m denotes the mass of final state W bosons. The Bl expressions for down type 
quarks are related to that of up type quarks as follows 

Bt{t,u,s) = Bl{u,t,s) , 



Bi{t,u,s) = B^{t,u,s) 



12) 



B^,{t,u,s) = -B^iu,t,s) . 



In addition, two more processes are allowed as the KK gravitons can appear at the 
propagator level, qq ^ G* WW and gg ^ G* ^ WW, as shown in Fig. US As 
we use unitary gauge in the electroweak sector the term proportional to 1/^ in the 
Vl^VT— graviton vertex |3] drops out. The qq and gg initiated contributions which are of 
order are given below. 



|M(o)|2 



1 



njsm^ - 16m^(t + u) + tu{'it^ + 2tu + 3^^) 



+m^(9t^ + mu + 9m^) - 2m^{t^ + It^u + Itu^ + u^)| 
-jsm® - 24m6(t + u) + tu{7f + lOtu + 7u^) 

+m^(17t2 + 62tu + 17^2) - 4m^{t^ + 9t\ + 9tu'^ + m^)| 



(13) 
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|M(o)|2 



99,9r 



(A^2 _ 1) 128 



X 



72 



-256m'*(t + u) + 192m^(t + - Q4m\t + uf '^^s^iAm^ - s) 

(n-l) 

, „ I - 44m*^ + 40m^(t + n) - AOm^tuit + u) + 9tu(t + uf 
[n — 2Y I 
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+m^(-9t^ + 26tu - 9m2) \ + ^^<| 692m^ - 13t^ - 196t^u - 362^^^ 

n — 2 



-imu^ - 13u^ - 544m^(t + u) - 8m'^{f + 83tu + 



u 



(14) 



We have denoted the sum of spin-2 KK graviton propagators by Vg, then times square 
of the couphng can be written as 



Stt 



d~2 



[rr + 2z/(A/v^)] 



(15) 



The function I{A/y/s) depends on the ultraviolet cutoff A on the KK modes and its 
expression can be found in [3]. The default choice for A would be the fundamental scale 
Ms unless mentioned otherwise. 

Next we give the interference of SM qq process with the gravity mediated qq subpro- 
cess. For convenience we will denote Mj^^^,,m^q-^]gr + c. by \M^'^^'^ ^- 



|M(o)|2 



e2 



qq,int 



N 



{cizi + cizi + cizi) 



(16) 



^0 



ryu 

^2 



V z rriz C^Tz cot dy/ XmVs 



[s — + T1 ml 



Zl = T^TZeVs 



TZeV, 



Cv Tz cot 9 



IS — mz 



w- 



+ 



2 Qu 



s — mlY + VI ml s 



(17) 



where are in their essence combinations of EW couphngs and propagator factors while 
Cl are the functions of kinematic invariants. These functions are given below 

= -^\Am^{-l + n)+t^it-u)u-m^{i-2Q + lln)t + nu} + 

mH {(-17 + 8n)t + (-11 + An) u} - m^t^^2 (-2 + n) f 

+ (-4 + n)tu + (-4 + n) 
1 



:i8) 



- lOm^ (-2 + n){t-u)+ (-17 + 8n) (t^ - u^) + 



8m2 

tn (t^ - u^) - 2m^ (-2 + n) [t^ - u^) \ , (19) 

Co" = -C^". (20) 
The Cf expressions for down type quarks are related to that of up type quarks as follows. 

Ct{t,u,s) = Cl{u,t,s) , 

Ci{t,u,s) = Cnt,u,s) , (21) 

c^{t,u,s) = c:;it,u,s) . 

Note that the W boson polarization sum —gfj,,y + k^j_ky/m? , which correctly takes into 
account 3-polarizations of a massive particle, does not give rise to negative powers of m 
and the m — )■ limit is smooth. 

2.2 Radiative Corrections 

In Fig. [ini the order loop diagrams that appear in SM and in Fig. [T2]the diagrams with 
a graviton propagator are presented. We will use Feynman gauge ^ = 1 in the QCD sector 
and we retain the term proportional to in the gluon-gluon-graviton vertex. However this 
term does not contribute to the matrix element squares. Here we consider only 5 flavors of 
quarks and treat them as massless. These diagrams contribute through their interference 
with the leading order diagrams. In general loop diagrams give ultraviolet divergences 
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and infrared divergences when the integration over loop momenta is carried out. We use 
dimensional regularization (n = 4 + e) to regulate these divergences; these divergences 
then appear as poles in e. Note however that owing to the gauge invariance and the fact 
that the KK gravitons couple to SM energy momentum tensor, a conserved quantity, 
this process is UV finite. The external parton leg corrections vanish in dimensional 
regularization for massless partons. From the loop Feynman diagrams in the appendix 
we find that all 2-,3-and 4-point loop integrals appear in the calculation. The maximum 
rank of tensor integrals is 3 and originate from the fermion box. These tensor integrals 
were reduced to scalar integrals following the procedure of Passarino- Veltman [21] . The 
4-point scalar integrals that appear in the gg initiated hox diagrams were taken from |25j . 
The one loop matrix elements are recorded below. The finite pieces of matrix element 
squares denoted by a superscript fin are given in the appendix. 
The SM contribution is found to be 



|M^|2 



qq,sm 



(^s{fiR)f{e,fi% s)Cb 



T(e) |M(o)|2 



qq,sm 



■fin 
qq,sm 



(22) 



the interference of SM with the gravity mediated processes are 



|M^| 



^qq,int = (^^j) /(c, /^L 



■fin 
qq,int 



IM^I 



gg,int 



■fin 
gg,int 



(23) 



(24) 



Note that in the above gg initiated SM diagrams with a quark triangle and 7/Z propa- 
gator do not contribute as they vanish in massless limit due to Furry' s theorem and weak 



isospin invariance. Also note that |M^|2^^ .^^ is completely finite; it does not contain any 
soft or collinear divergences because in SM the gg contribution begin at the loop level 
and a lowest order term should be finite. 
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The pure gravity contributions are 



q<l,9r 



99,9r 



as{lAi)f{^,IJ'R,s)CF 



Tf^(^)l^(°)|W + 4(2C(2)-5)|M(o)p 



qq,9r 



(25) 



16 4 



99,9r 



+ i (^72C(2) + 70^ - 203^ |M(o) |2 



(26) 



where 



(27) 



r(l + e) V4vr/x| 

The theory is renormahzed at scale /Xij. is the Casimir of the fundamental represen 
tation while Ca is the Casimir of adjoint representation in the color group. 



2N ' 



Ca = N, Tf = - 



(28) 



We can now write the order as{iJ,ji) contributions coming from virtual diagrams as, 
da"'-''^ = as{iiji)dxidx2f{e,ii%,s) 

--^ + — j ^ d(J^^l.{xi, X2, e) (/q,(xi)/g.(x2) + Xi O X2) 



^ , 16 4/11 
+Ca { + 



2 ■ CAe\^^^~ S^^^n J '^^SH^i^ ^2, e) (^fg{xi)fg{x2)^ 



+Cf ^ (i(jJ^^.*"(Xi, X2, e) {fqi{Xi)fq.{x2) + Xi ^ X2) 



+CAda^/'-{x„X2,e){f,{x,)f,{x2)) 



(29) 



The poles of order 2 in e in the one loop matrix elements correspond to the configurations 
which are both soft and coUinear simultaneously. These double poles cancel when real 
emission contributions are included, the remaining simple poles do not cancel completely 
and are factorized into the bare parton distribution functions at the scale /ip. 
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Several checks ensure the correctness of the matrix elements. The W boson polariza- 
tion sum —g^u + k^ku/m'^, which correctly takes into account 3-polarizations of a massive 
particle, does not give rise to negative powers of m and the m — )■ limit is smooth. Fur- 
ther, for gluon initiated process the gluon polarization sum is —g^u + {k^n^ + k^n^)/k.n 
where n is an arbitrary light like vector and the results are independent of the vector n. 
Furthermore the SM matrix elements are in agreement with the literature [8j. 

At NLO we also have to include 2 — )• 3 real emission processes. A generic process is 
of the form 

a{pi) + b{p2) ^ W^{p3) + W-{p,) + c{p,). (30) 

In Fig. [11] we show the qq and qg initiated real emission Feynman diagrams which appear 
in SM. In addition, in the ADD model the 2 — > 3 diagrams with graviton propagator are 
shown in Fig. [131 Here all the three kinds, qq,qg,gg initiated subprocesses occur. The 
2 — 7- 3 contributions to cross-section reveal the infrared divergences when the integral over 
the final state particles is carried out. The sum of virtual and real emission cross section 
is finite after mass factorization is carried out. For details we refer to the review [19] . 

Although the details of phase space slicing method to deal with soft and coUinear 
singularities in real emission processes were given in our earlier works p^[T^. we shall 
recapitulate it for completeness. The 2 — )• 3 phase space is divided into soft and coUinear 
regions using two small dimensionless slicing parameters 6s and 6c ■ The soft region 
is defined as the part of phase space where the final state gluon is soft and has an 
energy less than 6s^/s^2/'^ center of mass frame of incoming partons. The region 

complementary to the soft region is hard region and contains coUinear singularities. 
This region is thus further divided into hard coUinear region (the region of phase space 
where the final state parton is coUinear to one of the initial state parton) which contains 
coUinear singularities and hard non-colUnear region which is free of any singularities. 
All the order pieces together; the virtual cross-section da"""^ the soft piece da^°^^ and 
the mass f actor ized hard coUinear contribution da^^^^'^ (CT denotes mass factorization 
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counter term) is referred to as 2-body contribution. 

The only order piece, d(j^~^°'^y{5s-,5^^ which remains to be included is hard non 
coUinear and which is finite as the integration over 3-body phase space here does not 
include soft and collinear regions. The integration over the 3-body phase space is car- 
ried out using monte carlo, and it is constrained to avoid collinear and soft regions. 
The qq and gg initiated processes contain both kinds of divergences so the integral is 
constrained using 5s and 5c to avoid these regions. The qg initiated process, however, 
contain only collinear singularities (as soft fermions do not give singularities) and the 
3-body integration is constrained using only 5c. 

The NLO resuh is sum da^^ + da^-^°'^y{5s, 5c, HF) + da^-^"'^y{5s, 5c). The sum da^'^^'^y 
{5s, 5c, lip) +da^~'''°'^y{5s,5c) constitutes QCD correction, but da'^~^°'^y{5s,5c, fJ^p) and 
da^~''°'^y{5s,5c) independently are not physical quantities as these depend on the (arbi- 
trary) slicing parameters. The sum of these two pieces should be independent of the 
slicing parameters as these were introduced at the intermediate stages of calculation. 
We have checked that the sum of 2-body and 3-body contribution is independent of the 
5s and 5c over a large range of their values. 

In the next section we present the results using our monte carlo code which incorpo- 
rates the above given details. We will present the stability of results against variations 
of the slicing parameters. This code can easily accommodate any cuts on the final state 
bosons and can evaluate various kinematical distributions. 

3 Kinematical distributions and Results 

The LHC with a center of mass energy of 14 TeV will be our default choice. However 
we will also present some results for a center of mass energy of 10 TeV for the LHC. For 
numerical evaluation, the following SM parameters [26] are used 

mw = 80.398 GeV, mz = 91.1876 GeV, Tz = 2.4952 GeV, sin^ Ow = 0.231 (32) 
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where 6w is the weak mixing angle. For the electromagnetic coupling constant a we use 
= 128.89. CTEQ6 [27] density sets are used for parton distribution functions. 2-loop 
running for the strong coupling constant is used . The number of active massless-quark 
flavors is taken to be 5 and the value of Aqcd is chosen as prescribed by the CTEQ6 
density sets. At leading order, that is at order q;°, we use CTEQ6L1 density set (which 
uses the LO running a^) with the corresponding Aqcd = 165 MeV. At NLO we use 
CTEQ6M density set ( which uses 2-loop running ) with the Aqcd = 226 MeV; this 
value of Aqcd enters into the evaluation of the 2-loop strong coupling. The default choice 
for the renormalization and factorization scale is the identification to the invariant mass 
of the W boson pair ie., fip = f^R = Q- Furthermore the W bosons will be constrained 
to satisfy \yw\ < 2.5, where yw is the rapidity of a final state W boson . 
We will present below the following kinematical distributions: 

1. Invariant mass distribution, da/dQ, where Q is the invariant mass of the final state 
W boson pair, 

2. Rapidity distribution da/dY where Y = 1/2 ln(Pi ■ q)/{P2 ■ q), where Pi and P2 
are incoming proton momenta and q is the sum of the W boson 4-momenta. 

First we demonstrate that the sum of 2-body and 3-body contributions is fairly indepen- 
dent of the slicing parameters. In Fig. [1] (for SM) and Fig. [2] (for signal) we show the 
variations of these two pieces with the slicing parameters in invariant mass distribution 
at a value of invariant mass equal to 800 GeV. Here both 6s and 6c are varied together 
with the ratio 6s/ 6c fixed at a value of 100 [19]. We note that the sum of 2-body and 
3-body contributions is fairly stable against variations in these parameters and this gives 
us confidence in our code. In what follows we will use 6s = 10^^ and 6c = 10^^. 

In Fig. [3] we have plotted the invariant mass distribution both for the SM and the 
signal, in the range 300 GeV to 1300 GeV. In this plot we display for three extra 
dimensions ie., d = 3 and for fundamental scale equal to 2 TeV. To highlight the 
importance of QCD corrections we have also displayed the LO results of SM and the 
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signal, and we observe that the K factors (defined as = da^ / da^"^) are significantly 
large. We note that for the signal [sm + gr + int) K factor varies between 1.55 to 1.98 in 
the invariant mass range of 300 to 1300 GeV . This also shows that the LO results can 
be only treated as first approximations and to have more precise estimates we should 
go beyond the leading order. We note here that present computation does not take 
into account decay of W bosons to leptons which is observed experimentally, but as 
QCD corrections are independent of these decays, the K factors obtained here would 
not change when decays are taken into account. 

To estimate the effect of the number of extra dimension on the invariant mass distri- 
bution, we plot in Fig. Hlthe signal for three different values of d (3,4,5) with Ms fixed at 
2 TeV . We note that the lower the value of d, more is the strength of the signal. Next 
in Fig. Owe have plotted da/dQ for three different values of Ms (2.0, 2.5, 3.0) at a fixed 
value 3 for the number of extra dimensions. As expected, with increase in the funda- 
mental scale the deviations from SM predictions become less, and significant deviations 
from SM are observed at higher energies still. Next, in Fig. [6] we present the effect of 
variation of the UV cutoff A that appears in the expression of graviton propagator (see 
eq. [T5l) for Ms = 3TeV and d = 3. We observe that lowering the cutoff from Mg by 25% 
lowers the predictions by 14%. 

In Fig. [7] we have plotted the rapidity distribution da/dY at LO and NLO both for 
SM and the signal for d = 3 and Ms fixed at 2 TeV. We have plotted this distribution in 
the interval —2.0 < Y < 2.0 and have carried out an integration over the invariant mass 
interval 900 < Q < 1100 to increase the signal over the SM background. As expected 
the distribution is symmetric about Y = 0. 

As was noted before the NLO QCD corrections reduce the sensitivity of the cross 
sections to the factorization scale fip] this we now show in the Fig. |8l We have plotted 
SM and the signal both at LO and NLO, and have varied the factorization scale /ii? in the 
range Q/2 < fip < 2Q. The central curve in a given band (shown by the dotted curves) 
correspond to fip = Q- In all these the renormalization scale is fixed at fiR = Q. We 
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notice that the factorization scale uncertainties in SM are less compared to the signal. 
This is because of the dominant role of the gluon gluon initiated process in the signal. 
We see in this figure that a significant reduction in theoretical uncertainty, arising from 
the factorization scale, is achieved by our NLO computation. At Q = 1300 GeV the 
da/dQ for the signal varies by 18.8 % at LO as is varied between Q/2 to 2Q and it 
varies by 7.6 % at NLO. At the end we present in Fig. [9l da/dQ for LHC with a centre 
of mass energy of 10 TeV at NLO both for SM and signal. For comparison we have also 
plotted the 14 TeV results in the same figure. 

4 Conclusions 

In this paper we undertook computation of W boson pair production at the LHC at 
next-to-leading order in QCD in the extra dimension model of ADD. Here only spin-2 
KK gravitons appearing at the propagator level were considered. W boson production 
is one of the important channels at the LHC to probe both the standard model and 
new physics. As the leading order results serve only as first approximations we need to 
go beyond it to NLO to have more precise estimates. The NLO results are generally 
not only significantly larger as compared to the LO results but they are also much less 
sensitive to the arbitrary factorization scale and renormalization scale (if the LO already 
starts at order a^). 

Here we carried out a full NLO computation and presented analytical expressions 
of matrix element squares for all the SM, gravity mediated and the interference of SM 
and gravity mediated processes both at the LO and virtual level. We used dimensional 
regularization to regulate soft and coUinear divergences and the singularity structure in 
the loop level matrix elements is shown and it is observed that the singular pieces are 
proportional to the born contributions. As different kinematical distributions such as 
invariant mass distribution and rapidity distribution are evaluated with cuts on the final 
state W bosons it is useful to use monte carlo based semi analytical methods which 
allow to tailor code easily to these requirements. For this we used the two cutoff phase 
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space slicing method to divide phase space in soft and coUinear regions and filter out the 
singularities in the real emission contributions which appear on phase space integrations. 
A brief discussion on this method was presented. To save space we omitted the real 
emission matrix elements as the expressions are voluminous and can be obtained on 
request. We have used MS scheme throughout this paper. 

We have presented distributions for the LHC at lATeV and lOTeV. We first offered 
some checks on our monte carlo code such as stability of sum of 2 — body and 3 — body 
contributions against variation of the slicing parameters 6s and Sc and then presented 
invariant mass and rapidity distributions both at LO and NLO. We use CTEQ 6L1 and 
CTEQ 6M parton density sets for LO and NLO observables, respectively. Significant 
enhancements over the LO predictions are observed. The K factors are found to be 
large in the invariant mass distribution. For LHC at lATeV we find that K factor in 
the invariant mass distribution for the signal {sm + gr + int) varies between 1.55 and 
1.98 as Q varies between 300 to 1300 GeV. We have shown that a significant reduction 
in theoretical uncertainty, arising from the factorization scale, is achieved by our NLO 
computation. At Q = 1300 GeV the da/dQ varies by 18.8 % at LO as (Xp is varied 
between Q/2 to 2Q and it varies by 7.6 % at NLO. These observations justify the entire 
exercise and give results that are precise and suitable for further studies for constraining 
the parameters of the ADD model. Invariant mass distribution is also presented for LHC 
at a center of mass energy of lOTeV at the NLO level. 
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SM: 6g/6^=100 
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Figure 1: Variation of 2-body and 3-body contributions (of da/dQ at Q = 800 GeV 
SM) and their sum with 6s- Here 6s/ Sc = 100 has been used. 
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Signal: 83/8^=100 
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Figure 2: Variation of 2-body and 3-body contributions (of da/dQ at Q = SOOGeV^ in 
signal) and their sum with dg- Here 5s/ 5c = 100 has been used. 
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WW production at the LHC (V S =14 TeV) 
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WW production at the LHC (V S =14 TeV) 
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Figure 4: Effect of variation of number of extra dimensions in invariant mass distribution. 
The fundamental scale Mg has been fixed at 2 TeV. The curves correspond to NLO 
results. 
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WW production at the LHC (V S =14 TeV) 
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Figure 5: Effect of variation of tlie fundamental scale Mg in the invariant mass distri- 
bution. The number of extra dimensions has been fixed at 3. The curves correspond to 
NLO results. 
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WW production at the LHC (V S =14 TeV) 
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Figure 6: Effect of variation of tlie UV cutoff scale A in tlie invariant mass distribution 
for d = 3 and Ms = 3TeV. Tlie curves correspond to NLO results. 
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LHC (V S =14 TeV) 
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Figure 7: Rapidity distribution for Mg = 2TeV for SM and signal for d = 3 . We have 
integrated over the invariant mass range 900 < Q < 1100 to enhance the signaL 



23 



Hp variation, LHC (V S =14 TeV) 
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Figure 8: Factorization scale variation in the invariant mass distribution. The number 
of extra dimensions d = 3 and the fundamental scale Mg = 2TeV have been chosen. 
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LHC at V S =14 TeV and V S=10 TeV 
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Figure 9: Invariant mass distribution at NLO for SM and the signal. Here the thicker 
curves correspond to \/S = lOTeV and hghter curves to = UTeV at the LHC. 
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5 Appendix 



Below we give the finite pieces of the matrix element squares that appear in Eq. [22] 
through Eq. [2l]in section 2. 
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where contain the coupling and propagator factors , Bf are the functions of the 
kinematic invariants. 
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and the rest A^ are given in eqn. ([8]). The B^ functions are as follows 
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where contain the couphng and propagator factors while Cf are the functions of 
kinematic invariants. 

Zl = 'T^Tx XmV, (40) 
and rest are given in eqn. f|T7|) . The C? functions are as follows. 
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(i) (j) 

Figure 10: Leading order and order virtual diagrams in SM for the subprocess uu — )■ 
W^W" and gg — )■ W~^W~. The diagrams for the subprocess dd — )■ W~^W~ are obtained 
by replacing u d and W~ in the diagrams shown here. 
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(g) (h) 

Figure 11: Order real emission Feynman diagrams in SM for the subprocess uu — t- 
W^W~g and ug — )■ W^W~u. The diagrams for the subprocess dd — )■ W^W~g and 
dg — )■ W^W~d are obtained by replacing u ^ d and -H- iy~ in the diagrams shown 
here. 
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Figure 12: Leading order and order gravity mediated virtual correctionss for the 
subprocess qq — )■ W~^W~ and gg — )■ W^W~. The big dashed hues represent gravitons 
and the small dashed lines represent QCD ghosts. 
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Figure 13: Gravity mediated real emission diagrams for the subprocess qq — ?■ W^W g, 
qg — )■ W^W~u and gg — )■ W^W~g. The big dashed hnes represent gravitons. 
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